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We explore the phonon dynamics in thin films of the topological insulator material Bi 2 Se 2 Te using ultrafast 
pump-probe spectroscopy. The time resolved differential reflectivity in these films exhibit fast and slow 
oscillations. We have given a careful analysis of variation of phonon frequency as a function of film thickness, 
which we attribute to the existence of standing acoustic modes. However, no variation in the frequency of 
the optical phonon modes was found with film thickness. This indicates that the optical phonons intrinsically 
belong to the surface of topological insulators. The controllability of acoustic phonons by way of varying the 
film thickness will have tremendous implications in the application of these materials in low power spintronic 
device operating at room temperature. 
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The newly discovered 3D topological insulators (TIs) 
represent a class of exotic quantum state of matter with 
an insulating bulk and conducting surface states^ 
These conducting surface states exhibit a relativistic 
energy-momentum dispersion owing to the nontrivial 
Z 2 topology of the bulk electronic wave-function^. 
As a consequence, these surface states form a Dirac 
cone between bulk valence band and bulk conduction 
band which has been confirmed by angle resolved 
photo-emission spectroscopy (ARPES) studies^’®. Addi¬ 
tionally, the surface states are “topologically” protected 
from back-scattering. Helical spin texture due to high 
spin-orbit interaction, a tt Berry phase and the presence 
of odd number of Dirac cones are responsibe for this 
topological protection. Furthermore, the constraint 
of time reversal symmetry ensures that the surface 
electrons are not scattered from non-magnetic impurities 
and defect states. Apart from these, TI materials also 
exhibit interesting magneto-optical properties which 
arise due to the presence of axionic term®. These 
exciting properties have made TIs a subject of intense 
research both theoretically as well as experimentally 
over the past few years. TIs also promise important 
practical application in quantum computation and low 
power spintronic device applications^. 


The family of Bismuth based compounds including 
Selenides and Tellurides are well established examples 
of TI materials®. Previously, ultrafast pump-probe 
spectroscopy has been employed to study the phonon 
and carrier dynamics in Bi 2 Se 3 (BS)®’^®, Bi 2 Te 3 
(BT)^^’^^d3^ Sb 2 Te 3 ^^ etc. The differential reflectivity 
measured by pump-probe spectroscopy revealed fast 
and slow oscillations which have been attributed to 
surface optical and bulk acoustic phonons respectively^®. 
However these claims have not been conclusively sub¬ 
stantiated. In this letter, we reinvestigate these claims 


through similar single coloured, optical pump-probe 
spectroscopy in reflection geometry in TI films of differ¬ 
ent thicknesses. Most of the existing research have been 
carried out on BS and BT systems which have dominant 
bulk carriers owing to defect states that overwhelm 
the surface conduction^®. In our work, we have used 
Bi 2 Se 2 Te (BST) thin films, which exhibit lower bulk 
conduction owing to reduced defect states^^. Moreover, 
BST is also predicted to have a large spin texture 
giving rise to non-trivial Berry phase that prevents 
back-scattering^^’To support this claim, variation 
of the resistance with temperature has been shown in 
the supplementary material^®, clearly displaying the 
bulk insulating feature. We observed both fast and slow 
oscillations in BST thin films. 


BST thin films used in our experiments have been 
deposited on Silicon (100) substrates using pulsed laser 
deposition (PLD) technique. Samples of thickness 50 
nm, 100 nm and 300 nm were prepared in identical 
conditions by varying the number of laser pulses. Fur¬ 
ther details of sample preparation have been reported 
elsewhere^®. Furthermore, detailed information on the 
characterisation of these thin films are provided in the 
supplementary materiaP®. An ultrafast femtosecond 
pulsed laser (MANTIS from Coherent Inc. USA) was 
used to generate Gaussian pump and probe pulses of 
780 nm center wavelength and 100 fs pulse-width with 
a repetition rate of 80 MHz. The pump beam was 
modulated at 1.8 KHz using a mechanical chopper and 
pump induced change in the reflectivity of probe beam 
was measured using a lock-in amplifier as a function of 
pump-probe delay^®. The temporal resolution in our 
measurement was 66 fs (from 0 to 15 ps delay) and 333 
fs for the rest of the scan. Polarization of pump and 
probe pulse were made perpendicular to each other to 
reduce pump scattering. All measurements reported in 
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FIG. 1. Differential reflectivity of 100 nm thick BST film 
measured at room temperature is plotted (scattered plot) in 
semi-logarithmic scale. The solid line is the fitted curve. Inset 
shows the same data plotted in linear scale capturing the zero 
delay. 


this letter were performed at room temperature. 

Time resolved differential reflectivity in BST film is 
shown in fig.l (inset). A fast rise time indicates rapid 
transition of electrons from valence to conduction band. 
The semi-logarithmic plot of the data shows short 
lived rapid oscillations superimposed on rather slower 
oscillations that lasts for a longer duration. The fast 
oscillations last for around 10 ps, whereas, the slower 
one persists for around 100 ps. Measurements were 
carried for longer pump-probe delay (800 ps) but no 
significant modulation was observed after 100-120 ps. 

The differential reflectivity in fig.l is fitted with double 
decay function to capture the overall carrier relaxation 
dynamics. The damped oscillation terms were added to 
fitting equation to extract the characteristic frequencies 
along with the appropriate decay times^^. The fitting 
function is thus given by, 

Ai? 

—— = Ai exp {-(t - to)/ti} + A 2 exp {-(t - to)/t 2 } 

H 

+ Bi exp {-(t - to)As} sm(27r/i -fi (pi) 

+ B 2 exp{-(t - to)/t4} sin{2'Kf2 + (^ 2 ) + C (1) 

From the fit, the extracted frequencies of fast and 
slow oscillation are found to be 1.93 THz and 0.03 
THz respectively. These fast and slow oscillations are 
attributed to surface optical and bulk acoustic phonon 


FIG. 2. Raman scattering data (scattered plot) and its 
Lorentzian fit (solid line) of the BST thin film. Raman 
peak is at 65.6 cm“^ (1.97 THz). The other modes Eig and 
Afg are at 109.67 cm“^ and 158.74 cm“^ respectively. Inset 
shows Aig optical phonon mode vibration. Larger circles rep¬ 
resent Bi atoms wheras smaller circles represent Se and Te 
atoms. 


modes respectively^^. To confirm, Raman scattering 
spectra of these samples were taken using blue laser 
(488 nm) with a power of 6 mW and a magnification of 
50x. The charistic peaks are observed at 65.67 cm“^, 
109.67 cm“^ and 158.74 cm“^ (see fig. 2). Previous 
studies^^ attributed these peaks to the A}g, Eig and A\g 
Raman modes in BST (with hexagonal crystal structure) 
respectively. Frequency corresponding to the A\g Raman 
mode (1.97 THz) is in close agreement to the surface 
optical phonon mode observed in our pump probe 
measurements. The inset of fig. 2 shows the vibration 
corresponding to the A\g Raman mode. However, for 
the detection of lower symmetry Eg modes anisotropic 
polarization dependent measurements are necessary^^’^^. 
The Afg coherent phonons in Bismuth based systems 
have been reported to decay faster than A\g phonons 
and thus making them difficult to detect^^’^®. 


To understand the excitation density dependence, the 
pump fluence was varied from 15 /rJ/cm^ to 50 /rJ/cm^ 
while maintaining a constant probe fluence of 1.25 
/iJ/cm^. Frequencies of both the fast and slow oscilla¬ 
tions do not show any significant variation (see fig. 3). 
This indirectly supports the hypotheses that the two 
oscillations represent the phonon dynamics as change in 
the pump fluence will only change the carrier density 
in the conduction band and is not expected to have 
any effect on either of the phonon modes. However, we 
note that the amplitude of coherent oscillations changes 
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FIG. 3. Dependence of pump fiuence on fast (top panel) and 
slow oscillation (bottom panel) in BST thin film with a thick¬ 
ness of lOOnm. 



FIG. 4. Fitting (solid line) of differential reflectivity data 
(scattered plot) for BST samples of thickness (a) 50 nm, (b) 
100 nm, (c) 300 nm. 


almost linearly with fiuence. 


The same experiment was then repeated on identically 
prepared BST films of thicknesses 50 nm, 100 nm and 
300 nm. Any change in the thickness is expected to 
alter the boundary conditions for the acoustic phonon 
oscillations. In fig. 4, the pump-probe signal for the three 
different samples are shown along with their fitting. 
The magnitude of differential reflectivity is significantly 
reduced with the increasing thickness. To illustrate 
this further, the fast and the slow oscillations for all 
three samples are separately plotted in fig. 5 along with 
the ht. From this it can be easily concluded that both 


the frequency of oscillation and the damping factor for 
the fast oscillations remain fairly constant with varying 
thickness. However, this fast oscillation becomes more 
prominent in the thinner sample. Here, the damping 
factor of fast oscillations indicates the timescale of 
phonon phonon interactions^®. 

On the other hand, the slow oscillations are significantly 
affected by the change in thickness. The frequency of 
these oscillation decreases with increasing thickness. 
This can be explained by assuming that the slow 
oscillations are a result of interference between multiple 
reflection of the acoustic wave from the vacuum-sample 
and sample-substrate boundaries. The result of this 
interference is a standing wave like modulation which 
is expected to change with thickness as the boundary 
conditions are being effectively altered. This result is 
qualitatively similar to the normal mode oscillation of a 
string with fixed ends, where the frequency of oscillation 
is inversely proportional to the length of the string. 


Thickness of the sample also shows a pronounced 
effect on the damping of the slow oscillations. The 
damping can be physically attributed to the strain 
owing to the lattice mismatch between the sample and 
the substrate. Whereas hexagonal quintuple BST layers 
have a lattice constant of 4.218 and Si(lOO) has 
a diamond cubic lattice structure with lattice spacing 
of 5.431 A. Because of this lattice mismatch, there 
is a mechanical strain on the BST layers deposited. 
However, with increasing thickness, the BST films are 
allowed more degrees of freedom, which helps in relaxing 
this strain, resulting in reduced damping in the slow 
oscillations. From the fitting of fig. 4, the fast and and 
slow oscillations are extracted and plotted separately 
in fig. 5. The slow oscillations clearly demonstrate 
this reduction in damping with increasing thickness. 
Amplitude of the slow oscillations also decreases with 
increasing thickness. As mentioned above, the damping 
of the faster oscillations remains constant for different 
samples indicating that these are indeed characteristic 
of a surface phenomenon without any bulk contribution. 
Therefore, the optical phonons corresponding to the 
Ajg mode, intrinsically belong to the surface of the TI 
material. 

The ultrafast phenomena that induces optical phonon 
oscillations at the surface and acoustic phonon oscilla¬ 
tions in the bulk can be understood in the following 
way. The photo excited carriers experience a non-linear 
ponderomotive force at the surface^® and excite the 
cold lattice to oscillate coherently at the optical-phonon 
frequency (1.97 THz). The fast de-phasing of these 
oscillations are attributed to the existence of surface 
states^®. In the bulk, however, carriers relax through 
electron-acoustic phonon interaction process which is 
slower. Owing to the fact that the bulk is insulating 
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Sample 

Fast Oscillation 

Slow Oscillation 

Thickness 

Frequency 

Damping 

Frequency 

Damping 

(nm) 

(THz) 

(ps) 

(THz) 

(ps) 

50 

1.960 ± 0.008 

1.41 ± 0.03 

0.0330 ± 0.0008 

4.79 ± 0.21 

100 

1.914 ± 0.022 

1.40 ± 0.05 

0.0299 ± 0.0010 

8.42 ± 0.59 

300 

1.917 ± 0.049 

1.33 ± 0.03 

0.0149 ± 0.0010 

25.34 ± 1.03 


TABLE I. Observed frequency and damping of the fast and slow oscillations 


in these materials (BST) and the Fermi level resides in 
the bulk band gap, ultrafast dynamics of the carriers 
originating from the bulk is very different from that of 
surface^^, unlike in their parent compounds (BT and 
BS). 


From the perspective of device applications, though 
the surface electronic states are topologically protected 
against scattering from non-magnetic defects, acoustic 
phonons are still an impediment to the electron flow in 
the bulk. This is because, electron-phonon scattering 
will always be there at room temperature. For good 
device applications, in principle, one needs to lower the 
temperature to minimise the electron-phonon coupling. 
Moreover, the controllability of frequency of acoustic 
phonons by varying film thickness has a bearing on 
improving the low power spintronic device application 
and in certain cases, even without having to cool the 
system. 


In conclusion, the experiment captures the optical 
and acoustic phonon inodes through time resolved 
differential reflectivity measurements in BST thin films. 
Short lived fast oscillations, which are superimposed on 
the pump-probe signal are attributed to the optical 
phonon mode. The slower oscillations correspond to 
modulation due to interference of bulk acoustic phonon 
modes. By studying films of different thicknesses, we 
have been able to conclusively establish that optical 
phonon modes are indeed a surface phenomenon with 
no contribution form the bulk. On the other hand, 
the slow oscillations corresponding to the acoustic 
phonon modes are significantly affected by film thickness 
and this has been qualitatively explained in the light 
of dynamics of mechanical standing waves in elastic 
media. Our results will have significant impact on the 
application of topological insulators in spintronic device 
and topological quantum computation. 
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FIG. 5. Fast and slow oscillations for different samples as ex¬ 
tracted from fig.4. The dotted lines represent experimental 
data and solid lines represent the corresponding fit. (a)-(c); 
Fast oscillations in 50 nm, 100 nm and 300 nm films respec¬ 
tively. (d)-(f): Slow oscillations in 50 nm, 100 nm and 300 
nm films respectively. 
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